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■ Abstract 

OO 

We discuss the prospects for measuring the W mass in Run II of the Tevatron and at 
the LHC. The basic techniques used to measure Mw are described and the statistical, 
, . theoretical and detector-related uncertainties are discussed in detail. 

o 

m 

Measuring the W mass, Mw, is an important objective for the Tevatron experiments 
in Run II. One hopes to measure Mw with a precision of ±30 MeV per experiment per 
channel pQ. Together with the anticipated precision of 2 GeV for the top quark mass 
Ph. from Tevatron data [2., a determination of the W mass with an uncertainty of 30 MeV 
would make it possible to predict the Higgs boson mass, M H , from one-loop electroweak 
corrections to M w with an accuracy of approximately 30% (21 El Hj- With the best fit 
central value close to the LEP II direct search lower bound of Mh > 114.4 GeV |5J, and 
the rather poor \ 2 which results from a fit of the current precision electroweak data to 
the Standard Model (SM) [3IZj, this would bring considerable additional pressure on the 
SM. 

This document is structured as follows. The basic techniques used to measure the 
W mass and the expected statistical and overall uncertainties are briefly discussed in 
Sec. El The detector-related and theoretical uncertainties affecting the W mass and width 
measurements are described in more detail in Sec. H3 and Sec. 01 respectively. Alternative 
methods of measuring the W mass at the Tevatron are discussed in section El Section El 
concludes this document. 

2 M w Measurement from the M T Lineshape 

The determination of Mw depends on the two body nature of the W decay, W — * ivg. 
The kinematical Jacobian peak and sharp edge at the value of Mw /2 is easily observed in 

1 Invited talk presented at the Mini- Workshop "Electroweak Precision Data and the Higgs Mass", 
DESY Zeuthen, Germany, February 28 - March 1, 2003, to appear in the Proceedings 
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Figure 1: The effects of resolution and 
the finite on pr(e) in W boson decay. 
The histogram shows p^ without detector 
smearing and for pjf = 0. The dots in- 
clude the effects of adding finite Pr , while 
the shaded histogram includes the effects 
of detector resolutions. The effects are cal- 
culated for the D0 Run I detector resolu- 
tions. 
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Figure 2: The effects of resolution and the 
finite Pr on Mt in W — > ev. The his- 
togram shows Mt without detector smear- 
ing and for pYf = 0. The dots include 
the effects of adding finite pjf , while the 
shaded histogram includes the effects of 
detector resolutions. The effects are cal- 
culated for the D0 Run I detector resolu- 
tions. 



the measurement of the transverse momentum (jpr) of either of the leptons. However, the 
Jacobian peak in the px distribution is smeared out by finite W width effects, detector 
resolution effects, and, in particular, the effect of the finite transverse momentum of the 
W boson, pYf . Figure H] shows a calculation of the Px( e ) distribution (unsmeared) with 
p^ = 0; the effect of finite p^; and the inclusion of detector smearing effects. It is 
apparent that pr(e) is very sensitive to the transverse motion of the W boson. 

Historically, precise understanding of p^ has been lacking, although it is currently 
modeled by measurable parameters through the resummation formalism of Collins, Soper, 
and Sterman |8 a . For this reason, the transverse mass, Mt, was suggested [9^ and has 
been the traditional measurable. It is defined by 

M T = ^/2p T (£K(z/)(l-cos(0^)), (2.1) 

where <fit tV is the angle between the charged lepton and the neutrino in the transverse 
plane. The observables are the lepton transverse energy or momentum pr(^) and the 
non- lepton transverse energy u (recoil transverse energy against the W), from which the 
neutrino momentum pr( z/ ) and the transverse mass My are derived. Mt is invariant 
under transverse boosts to first order in the velocity of the W boson. Figure El shows that 
the sensitivity of Mt to p^ is indeed almost negligible. While considerably more stable to 
the phenomenology of the production model, the requirement that the neutrino direction 
be accurately measured leads to a set of experimental requirements which are difficult in 
practice to control. 

Both CDF and D0 have determined the W boson mass using the transverse mass 
approach. The combined pp collider result is [7j 

M w = 80.454 ± 0.059 GeV. (2.2) 
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Figure 3: Statistical uncertainties in Run I Mw measurements. Each circle represents 
either a CDF or D0 measurement. The result of a straight line fit is shown. The shaded 
box is the approximate extrapolation to a 2 fb -1 Run II result. 

Presuming that Run II is to deliver an integrated luminosity of 2 fb _1 , the statistical 
precision on Mw can be estimated from the existing data. Figure |3] shows the statis- 
tical uncertainties in these measurements as a function of 1/y/Nw (Nw is the number 
of W boson events), demonstrating a predictable extrapolation to Nw — 700,000 which 
corresponds to the expected Run II dataset per experiment per channel. The statisti- 
cal uncertainty from this extrapolation is approximately 13 MeV. Taking into account 
the detector related and theoretical uncertainties discussed in more detail below, a total 
uncertainty of ±30 MeV for Mw is projected. The current CDF and D0 Run II data 
samples are approximately 70 pb _1 [TU] . 

For integrated luminosities larger than 2 fb -1 , the systematic errors rapidly become 
dominant. For 15 fb _1 , the goal for Run lib, one hopes to achieve a precision of [TT] 

5M W = ±17 MeV. (2.3) 

For even larger data sets, the determination of the energy resolution will be systematically 
limited by the uncertainty in the width of the Z boson [llj. This will restrict the accuracy 
to about 15 MeV per channel per experiment. 

At the LHC, for 10 fb _1 , one expects about 6 x 10 7 W — ► Iv (£ = e, /x) events, resulting 
in a statistical error of < 2 MeV • If the lepton energy and momentum scales can be 
determined with a precision of 0.02% or better, the W mass can be measured with a total 
uncertainty of ±15 MeV H2J - 
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Figure 4: Systematic uncertainties for each Run I My/ measurement. The open squares 
are the four electron measurements from CDF and D0, the circles are the scale uncer- 
tainties from two D0 electron measurements and the Run lb CDF measurement, and the 
diamonds are the systematic uncertainties (excluding scale) for the CDF muon measure- 
ments. The large box is the position of the extrapolated statistical uncertainties to the 
Run II luminosity. The lines are linear fits to each set of points. 

3 Detector-specific Uncertainties 

Much of the understanding of experimental systematics comes from a detailed study 
of Z bosons and hence as luminosity improves, systematic uncertainties should diminish 
in kind. Certainly, the scale and resolution of the recoil energy against the W come 
from measurements of the Z system. Likewise, background determination, underlying 
event studies, and selection biases depend critically, but not exclusively, on Z boson 
data. Most importantly, the lepton energy and momentum scales depend solely on the Z 
boson datasets. In the Run I W mass measurement, the largest systematic uncertainties 
originated from the lepton energy and momentum scales, and the modeling of the W 
recoil. Non-Z related recoil systematics were estimated to enter at the 10 MeV level. 

Figure 0] shows the CDF and D0 systematic uncertainties for both electrons and 
muons as a function of 1/y/Nw- The calorimeter scale uncertainties for electrons is cur- 
rently tied to the determination of fiducial di-lepton decay resonances, notably the Z 
boson, but also the J/ip, T and the E/p dependence on the energy E, using electrons 
from W and Z decays. As the statistical precision improves, the dominant feature of 
the scale determination becomes its value in the region of Myy, so offsets and any low 
energy nonlinearities become relatively less important and hence reliance on the low mass 
resonances is reduced. On the other hand, for the muon momentum scale determination, 
where the observable is the curvature, low mass resonances are also important. Figure 0] 
suggests that this uncertainty is truly statistical in nature and extrapolates to approxi- 
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mately the 15 MeV level. The ability to bound non-linearities using collider data may 
become a limiting source of error in Run II. 

Figure 0] also shows the non-scale systematic uncertainties from both the CDF and 
D0 electron measurements of M\y and the CDF muon measurement. Here the extrapo- 
lation is not as straightforward, but there is clearly a distinct statistical nature to these 
errors. That they appear to extrapolate to negative values suggests that the systematic 
uncertainties may contain a statistically independent component for both the muon and 
the electron analyses. 

Z — > £ + £~ data constrain both the lepton scales and resolutions and an empirical model 
of the hadronic recoil measurement. The determination of the lepton scales and resolutions 
requires a measurement of the Z boson mass, Mz and width, Tz, and calibration of these 
observables using the available LEP data. In order to do this, it is necessary to use 
exactly the same theoretical input which has been used to extract Mz and Tz at LEP, ie. 
to include the reducible and irreducible 0(g 4 mf jMy V ) corrections to the effective leptonic 
weak mixing parameter, sin 2 ^, and the W mass ^3], hi the calculation. A calculation 
of the complete 0(a) corrections to pp — > 7, Z — > £ + £~ which also takes into account 
the 0{g A m^/M^ v ) corrections has recently been completed [Tlj . 

4 Theoretical Uncertainties 

The Mt lineshape simulation requires a theoretical model, as a function of the W 
mass and width, of the differential cross section. The W mass is extracted by comparing 
data with templates of the Mt distribution for different values of Mw- The theoretical 
model must take into account both QCD corrections and electroweak corrections. 

QCD provides an understanding of the p^ distribution and thus of the hadronic recoil 
in W events. For a quantitative understanding of the W (and Z) boson pt distribution 
it is necessary to resum soft gluon radiation. In the low pt region, the perturbative 
calculation must be augmented by a non-perturbative contribution which depends on 
three parameters which are tuned to fit the Z — > £ + £~ data. Electroweak radiative 
corrections, in particular final-state QED radiation is important, because bremsstrahlung 
shifts the W mass to smaller values and affects the isolation variables needed to select a 
clean W sample. 

To date, ad hoc event generators have been used in the W mass measurement. In 
Run II, this measurement will reach a precision of tens of MeV, requiring much more 
attention to detail in Monte Carlo calculations. Precision electroweak measurements in 
Run II should strive to use Monte Carlo programs that are common to both collider 
experiments and include state-of-the-art QCD and electroweak calculations. Currently, 
simulations are largely based on two programs: RESBOS jT3] which includes resummed 
QCD corrections but ignores electroweak radiative corrections, and WGRAD jTH] which 
includes 0(a) electroweak corrections, but ignores QCD corrections. A merging of the 
two programs would be highly desirable in order to achieve the goal of a unified program 
which includes both QCD and electroweak radiative corrections. 

The current W production model has a number of shortcomings which lead to sys- 
tematic uncertainties in the W mass measurement. These uncertainties are discussed in 
the following. 
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4.1 Parton Distribution Functions 



The transverse mass distribution is invariant under the longitudinal boost of the W 
boson. However, the incomplete rapidity coverage of the detectors introduces a depen- 
dence of the measured Mr distribution on the longitudinal momentum distribution of 
the produced W's, determined by the parton distribution functions (PDF's). Therefore, 
quantifying the uncertainties in the PDFs and the resulting uncertainties in the W mass 
measurement is crucial. 

The measurement of the W charge asymmetry at the Tevatron, which is sensitive to 
the ratio of d to u quark densities in the proton, is of direct benefit in constraining PDF 
effects in the W mass measurement. This has been demonstrated by the CDF experiment. 
From the variation among the six reference PDFs, an uncertainty of 15 MeV was taken 
which is common to the electron and muon analyses. Since the Run lb charge asymmetry 
data is dominated by statistical uncertainties, a smaller uncertainty is expected for the 
Run II measurement. Measurements of Drell-Yan production at the Tevatron can be used 
to get further constraints on PDFs. 

Since the PDF uncertainty originates from the finite rapidity coverage of the detectors, 
it is expected to decrease with the more complete rapidity coverage of the Run II detectors. 
The advantage of a larger rapidity coverage has been demonstrated by the D0 experi- 
ment ^7j: the uncertainty on the W mass measurement using their central calorimeter 
was 11 MeV, while that using both the central and end calorimeter was 7 MeV. With the 
upgraded calorimeters and trackers for the rapidity range \r}\ > 1, the CDF experiment 
can measure the W mass over a larger rapidity range in Run II. 

There has been a systematic effort to map out the uncertainties allowed by available 
experimental constraints, both on the PDFs themselves and on physical observables de- 
rived from them ^HJEI]- This approach will provide a more reliable estimate and may 
be the best course of action for precision measurements such as the W mass or the W 
production cross section. 

4.2 W Boson Transverse Momentum 

The neutrino transverse momentum is estimated by combining the measured lepton 
transverse momentum and the W recoil: Pt{v) = — {pt{C-) + u). It is clear therefore that 
an understanding of both the underlying W boson transverse momentum distribution 
and the corresponding detector response, usually called the recoil model, is crucial for 
a precision W mass measurement. For the CDF Run lb W mass measurement, the 
systematic uncertainties from these two sources were estimated to be 15 — 20 MeV and 
35 — 40 MeV, respectively, in each channel |20j . 

The strategy employed in Run I, which is expected to be used also in Run II, is to 
extract the underlying pYji distribution from the measured distribution (Y is the weak 
boson rapidity): 

d 2 a d 2 a ^ d 2 a/dpfdY 

dp^dY ~ dp^dY X d?a/dp^dY ' ^ ' ' 

where the ratio of the W and Z differential distributions is obtained from theory. This 
method relies on the fact that the observed distribution suffers relatively little from 
detector smearing effects, allowing fits to be performed for the true distribution. The 
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Figure 5: (a) A comparison of the two sources of uncertainty on the derived pj. distribu- 
tion, (b) The p^ distribution extracted for the CDF Run lb W — > \xv mass measurement. 

experimental uncertainties, as in many aspects of the W mass measurement, are dom- 
inated by the available Z — > £ + £~ statistics and should scale correspondingly with the 
delivered luminosity in Run II. Theoretical uncertainties in the ratio ofWtoZ transverse 
momentum distributions contribute a further 0(5) MeV to the overall error. The two 
sources of uncertainty are compared for the CDF Run lb W — > \w analysis in Fig. 03 

4.3 QCD Higher Order Effects 

The charged leptons in W — > iv are produced with an angular distribution determined 
by the 0(a 2 s ) calculation of Ref. [21] which, for W + bosons with a helicity of —1 with 
respect to the proton direction, has the form: 

, oc 1 + a x ($) cosfl C s + a 2 (pW) cos 2 fl cs . (4.2) 

CL COS Uq§ 

Here, 8qs is the polar direction of the charged lepton with respect to the proton direction 
in the Collins-Soper frame [22] . a± and a 2 are pt dependent parameters. For p^ = 0, 
ai = 2 and a 2 = 1, providing the angular distribution of a W boson fully polarized along 
the proton direction. For the p^ values relevant to the W mass analysis {p^ < 30 GeV), 
the change in W polarization as p^ increases only causes a modest change in the angular 
distribution of the decay leptons [2*T] . 

While the uncertainty associated with the change in the angular distribution of the 
W decay lepton due to higher order QCD corrections (a few MeV) has been negligible for 
the Run I measurements, it can not be ignored for the Run II measurements. 



7 



4.4 Electroweak Radiative Corrections 

The understanding of electroweak radiative corrections is crucial for precision W mass 
measurements at the Tevatron. The dominant process is final state photon radiation 
(FSR) from the charged lepton, the effect of which strongly depends on the lepton iden- 
tification criteria and the energy or momentum measurement methods employed. Calori- 
metric energy measurements, such as those employed in the electron channel, are more 
inclusive than track based momentum measurements used in the muon channel and the 
effect of FSR is consequently reduced. In the CDF Run lb W mass measurement the mass 
shifts due to radiative effects were estimated to be —65 ±20 MeV and —168 ±10 MeV for 
the electron and muon channels, respectively ,20j. These effects will be larger in Run II 
due to increase in tracker material in CDF and magnetic tracking in D0. 

The Monte Carlo program used for the Run I W mass measurement incorporated 
a calculation of QED corrections by Berends and Kleiss This treatment, however, 
does not include initial state radiation (ISR), or soft and virtual corrections. Furthermore, 
multi-photon radiation effects are ignored in Ref. j2H|. The systematic uncertainties of the 
missing radiative corrections were estimated to be 20 MeV and 10 MeV in the electron 
and muon channels, respectively [2U] ■ Clearly these systematic uncertainties become much 
more significant in the context of statistical uncertainties of (9(10) MeV expected for 2 fb _1 
in Run II. 

In the last few years two independent calculations of the full 0(a) electroweak (EW) 
corrections to W production in Monte Carlo form have been carried out ^Hl I2U 121] • Cal- 
culating the EW radiative corrections to resonant W boson production, the problem arises 
how an unstable charged gauge boson can be treated consistently in the framework of per- 
turbation theory. This problem has been studied in Ref. [26J with particular emphasis on 
finding a gauge invariant decomposition of the EW 0(a) corrections into a QED-like and 
a modified weak part. In Ref. [21] 5 it was demonstrated how gauge invariant contributions 
that contain the infrared (IR) singular terms can be extracted from the virtual photonic 
corrections. These contributions can be combined with the also IR-singular real photon 
corrections in the soft photon region to form IR-finite gauge invariant QED-like contribu- 
tions corresponding to initial state, final state and interference corrections. The collinear 
singularities associated with initial state photon radiation can be removed by universal 
collinear counter terms generated by "renormalizing" the parton distribution functions, 
in complete analogy to gluon emission in QCD. 

The various individual contributions to the electroweak 0(a) corrections on the Mt 
distribution are shown in Fig. |U1 The initial state QED-like contribution uniformly in- 
creases the cross section by about 1% for electron (Fig. HJi) and muon (Fig. EJd) final 
states. It is largely canceled by the modified weak initial state contribution. The in- 
terference contribution is very small. The final state QED-like contribution significantly 
changes the shape of the transverse mass distribution and reaches its maximum effect in 
the region of the Jacobian peak, Mt ~ Mw- 

Since final state photon radiation has a significant impact on the W mass extracted 
from data, one has to worry about the effects of multi-photon radiation. A first attempt 
towards the inclusion of the 0(a 2 ) real photon corrections in qq' — > W — » iwyy and 
qq — > Z — > ^ + £~77 was made in Ref. More recently, two calculations of higher 
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Figure 6: Ratio of the 0(a 3 ) and lowest order cross sections as a function of the transverse 
mass for a) pp — > e + u{^i) and b) pp — > fj, + v(j) at y/s = 1.8 TeV for various individual 
contributions. The upper (lower) solid lines show the result for the QED-like initial (final) 
state corrections. The upper (lower) dotted lines give the cross section ratios if both the 
QED-like and modified weak initial (final) state corrections are included. The dashed lines 
display the result if only the initial - final state interference contributions are included. 
Lepton identification requirements are not taken into account here. 

order QED corrections to W — > Iv decays have been performed which include both real 
and virtual corrections. Ref. [2H] uses Yennie-Frautschi-Suura exclusive exponentiation 
to calculate multi-photon radiation effects whereas a QED parton shower approach is 
employed in Ref. j2H]- Using a simplified detector model, multi-photon radiation in W 
decays was found j2S| to shift the W mass by a few MeV in the electron channel, and by 
about 10 — 15 MeV in the muon channel. In view of the expected precision of 15 — 30 MeV 
for Mjy in Run II and at the LHC, it will be important to take multi-photon radiation 
effects into account when extracting the W mass from data. 

So far, no calculation exists which includes multi-photon radiation effects in hadronic 
Z boson production. In Z — » both final state leptons couple to photons, in contrast 

to W — > iv where only one of the final state particles carries electric charge. As a result, 
one expects that multi-photon radiation affects the Z boson mass and width more strongly 
than the corresponding W observables. Since Mz and Y z are crucial in determining the 
lepton energy and momentum scales, a calculation of Z boson production in hadronic 
collisions which takes into account higher order QED corrections is needed. 
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5 Other Methods of Determining M w at Hadron Colliders 

While the traditional transverse mass determination has been the optimal technique 
for the extraction of My/ in previous hadron collider experiments, other techniques may 
be employed in the future. These methods may shuffle or cancel some of the systematic 
and statistical uncertainties resulting in more precise measurements. 

5.1 Transverse Momentum Fitting 

As noted in Sec. 121 the most obvious extensions of the traditional transverse mass 
approach to determining My/ are fits of the Jacobian kinematical edge from the transverse 
momentum of both leptons. D0 has measured My/ using all three distributions and the 
uncertainties are indeed ordered as one would expect: The fractional uncertainties on My/ 
from the D0 Run I measurements for the three methods of fitting are: 0.12% (Mr), 0.15% 
(pr(e)), and 0.21% (p T (i/)). As expected, the pr(e) method is slightly less precise than 
the measurement using the transverse mass. However, for a central electron (|^| ;$ 1), 
the uncertainty in the pr(e) measurement due to the Pr model is 5 times that in the 
Mt measurement. This is nearly balanced by effects from electron and hadron response 
and resolutions which are relatively worse for Mt- Accordingly, when there are sufficient 
statistics to enable cuts on the measured hadronic recoil, the measurement uncertainty 
from the model might be better controlled and enable the pr(e) measurement to 
compete favorably with the Mt measurement which relies so heavily on modeling of the 
hadronic recoil. In order to optimize the advantages of all three measurements, the D0 
final Run I determination of My/ combined the separate results [30] . 

The resolution sensitivity for muon measurements is even less than that for electrons 
so that has the benefit of slightly favoring a transverse mass measurement with muons 
over that for electrons. 



5.2 Ratio Method 

The W mass can also be determined from the ratio of the transverse mass distributions 
of the W and Z boson |3TJ E2 EH] • The advantage of the method is that one can cancel 
common scale factors in ratios and directly determine the quantity r meas = which 
can be compared with the precise value of Mz from the LEP experiments. The downside 
of the ratio method is that the statistical precision of the Z sample is directly propagated 
into the resultant overall 5 My/- 

D0 has carried out a "proof of principle" analysis, determining My/ from the ratio of 
the transverse mass distributions of the W and Z boson jHH E21 EH] • Since the W and 
Z boson mass are different, the Z boson data are rescaled to fit the W boson data. The 
preliminary D0 result for My/ from Run lb data using the transverse mass ratio method 

is (23 

M w = 80.115 ± 0.211 (stat.) ± 0.050 (syst.) GeV. (5.1) 
For comparison, the result from the fit to the Mt distribution is jHO] 

M w = 80.440 ± 0.070 (stat.) ± 0.096 (syst.) GeV. (5.2) 
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Table 1: Projected uncertainties of the W boson mass measurement using the transverse 
mass fit and the transverse mass ratio method for W — > ev. 



J Cdt (fb- 1 ) 2 15 30 


/ Cdt (fb- 1 ) 2 15 30 


M T fit 

stat. uncert. (MeV) 19 7 5 
syst. uncert. (MeV) 19 16 15 


ratio method 

stat. uncert. (MeV) 44 16 11 
syst. uncert. (MeV) 10 4 3 


total uncert. (MeV) 27 17 16 


total uncert. (MeV) 45 16 11 



As expected, the transverse mass ratio method results in a much improved systematic 
uncertainty. The statistical uncertainty, however, is significantly larger than for the fit to 
the Mr distribution. 

As mentioned in Sec. |2l systematic uncertainties begin to limit the precision for M\y 
which can be achieved from the fit to the transverse mass distribution at the Tevatron for 
integrated luminosities of 15 fb" 1 or more. Such high luminosities may be achievable in 
Run lib. Up to 30 fb -1 are conceivable if additional upgrades are made to the accelerator 
complex [TT]. Table ^ shows the expected precision of the W mass measurement from 
the transverse mass fit and the W to Z transverse mass ratio, extrapolated from the 
Run lb measurement by D0 jHOj- For integrated luminosities larger than about 15 fb -1 
the transverse mass ratio method gives an uncertainty for Mw which is smaller than that 
obtained with the traditional method of fitting the distribution. 

6 Conclusions 

The measurements of the W mass in Run I of the Tevatron contribute significantly to 
the world average of Mw- Close inspection of the various systematic error sources leads 
us to believe that a W mass measurement in Run II at the 30 MeV level per experiment is 
achievable. Alternative methods for determining Mw at the Tevatron have been discussed 
and may turn out to be more appropriate in the Run lib operating environment than the 
traditional transverse mass fitting approach. Determination of the W mass at the LHC 
will be extremely challenging, using detectors that are not optimized for this measurement. 
Clearly, the W mass measurement at the Tevatron in Run II will contribute significantly to 
a reduction of the overall uncertainty of Mw, and thus to a precise indirect determination 
of the Higgs boson mass. 
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